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In this work we analyze a non-interacting one dimensional polymer Bose-Einstein condensate in 
an harmonic trap within the semiclassical approximation. We use an effective Hamiltonian coming 
from the polymer quantization that arises in loop quantum gravity. We calculate the number of 
particles in order to obtain the critical temperature. The Bose-Einstein functions are replaced 
by series, whose high order terms are related to powers of the polymer length. It is shown that 
the condensation temperature presents a shift respect to the standard case, for small values of the 
polymer scale. In typical experimental conditions, it is possible to establish a bound for A 2 up to 
< 10~ 16 m 2 . To improve this bound we should decrease the frequency of the trap and also decrease 
the number of particles. 

PACS numbers: 04.60Bc, 04.60.Kz, 04.60.Pp, 03.75.Nt 



I. INTRODUCTION 

Several quantum gravity models suggest that the rela- 
tion between energy and momentum of microscopic par- 
ticles must be modified as a consequence of the quantum 
structure of space-time [J-Q • A deformed dispersion re- 
lation, emerges as an ideal tool in the search for possible 
effects related to the quantum structure of space-time. 
Nevertheless, the most difficult aspect in the search of 
experimental evidence relevant for the quantum gravity 
problem is the smallness of the possible effects [H |j| • Un- 
fortunately, due to this fact, the possible bounds for the 
deformation parameters, open a wide range of possible 
magnitudes, which implies a significant challenge. 

Among the models that introduce such deformations, 
polymer quantum systems are simple quantum mechan- 
ical models quantized in a similar way as in loop quan- 
tum gravity [H [H|. These systems provide scenarios 
where some characteristics of the full theory, specially 
the discrete nature of space, can be explored in a simpler 
context [l2l-[l6j]. Indeed, in polymer quantum mechan- 
ics the momentum operator can not be defined, hence a 
regularized operator is proposed b y t he introduction of 
the so-called polymer length scale [1(| [H| • Particularly, 
from this model an effective Hamiltonian which contains 
some trace of the granularity of space can be drawn [llj ■ 

The use of Bose-Einstein condensates as a tool in the 
search of quantum gravity effects, for instance, in the 
context of Lorentz violations, or to provide phenomeno- 
logical constrains on Planck scale physics, has produced 
several interesting research [j| [(| [l8l425j (and references 
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therein). In these works, the possible effects arising from 
Planck scale physics by looking at some modifications 
in the thermodynamic properties of Bose-Einstein con- 
densates was explored. Since its observation with the 
help of magnetic traps (27rl30| , the phenomenon of Bose- 
Einstein condensation, from the experimental and theo- 
retical point of view, has spurred an enormous amount of 
publications (30M5TI l58l. 159] (and references therein). In 
particular the condensates have been studied in different 
spatial dimensions [13, [H2, [H, [59| . Although, a one di- 
mensional condensate can never be reached, it is possible 
to obtain a quasi one dimensional condensate, by using 
extremely anisotropic traps [13, [H, H3, [H| . Additionally, 
the condensation process of a bosonic gas in one dimen- 
sion, trapped in a box or in an harmonic oscillator po- 
tential, has a very particular behavior. In this situation, 
apparently the condensation is not possible in the ther- 
modynamic limit if one assumes that the ground state 
energy is zero, or equivalently, if the associated chemical 
potential is zero at the condensation temperature. How- 
ever, taken into account the minimum energy associated 
to the system, the condensation process is possible at 
finite temperature |H EI H • 

The aim of this work is to analyze the properties of 
a polymer one dimensional condensate trapped in a har- 
monic oscillator potential, within the semiclassical ap- 
proximation, in order to obtain representative bounds 
for the polymer length scale. We will show also that for 
this one dimensional system the functional form of the 
relevant thermodynamic quantities, particularly, for the 
number of particles, can be expressed by series of the 
so-called Bose-Einstein functions. This series converges 
when the polymer scale is small enough. Finally, by the 
introduction of the ground state energy associated to this 
system, we calculate the shift in the usual condensation 
temperature caused by the polymer scale A 2 . These facts 
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allow us to bound the polymer scale by using one dimen- 
sional finite size systems. 



II. POLYMER QUANTIZATION 

Let us now present some of the main results of polymer 
quantization [10 . 1 111 . This quantization arises from loop 
quantum gravity [561 ] . In the loop or polymer representa- 
tion the corresponding Hilbert space "Hpoiy is spanned by 
the basis states {|xj)}, whose coefficients have a suitable 
fall-off [lfj, with the following inner product 



(xi\xj) = 5, 



(1) 



where S^j is the Kronecker delta. The polymer Hilbert 
space can be represented as H po i y = L 2 (Rd,d[id), where 
d\id is the corresponding Haar measure, being R<j the real 
line endowed with the discrete topology, that is, the dual 
of the Bohr compactification of the real line [lfl, [56[ . The 
basic operators in this quantization are the position and 
translation. The position operator x acts as usual by 
multiplication 

while the translation operator V(X) moves to a position 
of arbitrary distance A 



V(\)\ Xj 



A). 



In the Schrodinger quantization, the operator V(A) is 
weakly continuous in A and the momentum operator is 
defined as its infinitesimal generator. However, in the 
polymer representation the translation operator V(A) 
fails to achieve this condition, i.e., it is not weakly con- 
tinuous on A due to the discrete structure assigned to 
space p.0]. Therefore, there is no Hermitian operator as 
infinitesimal generator of translations and thus, the mo- 
mentum operator is not well defined. Due to this fact, 
any phase space function has to be regularized by the 
introduction of a lattice of fixed positive length. We in- 
terpret this length as a fundamental scale that is called 
the polymer length scale [lfj. With this regularization, 
the polymer Hamiltonian can be defined as follows 



H x = 



2mA 2 



V(X)-V(-X) +U(x) 



(2) 



where U(x) is the potential term. The action of the 
Hamiltonian ([2} decomposes the polymer Hilbert space 
Wpoiyj m t° a continuum of separable superselected sub- 
spaces, each with support on a regular lattice {nA + 
Xo\n £ Z}, with xq € [0, A) that parameterizes a particu- 
lar superselected sector [1J| ■ Notice that the Hamiltonian 
([2]) can be formally written as 



2mA 2 ' 




U(x). 



(3) 



We can use this expression to obtain the effective Hamil- 
tonian, simply by replacing the kinetic term with the 
square of the sine function. 

It is worth mentioning, as was pointed out in [17| , that 
the effective three dimensional case present some compli- 
cations in calculating analytically the integrals. Due to 
this fact, in the present work we restrict ourselves to the 
one-dimensional case. 



III. CONDENSATION TEMPERATURE 

Let us start with a one-dimensional effective polymer 
Hamiltonian 17[, given by 



H = 



2mA 2 bm V h 



U(x) 



(4) 



where U(x) = muj^x 2 /2 is the one dimensional harmonic 
oscillator potential that model the trap, and A is the so- 
called polymer length scale. Therefore, the semiclassical 
energy spectrum associated to ^ can be expressed as 
follows 



(t) 



2mA 2 



(5) 



The associated one-dimensional spatial density can be 
written as [34|, l59j 



n(x) 



I 

2Trh 



dp x 



o/3(e p -M) _ 1 ' 



(6) 



where, as usual, (3 — 1/kT, being k the Boltzmann con- 
stant, and /i is the corresponding chemical potential. 
Consequently, the total number of particles of the sys- 
tem is 



N = n(x)dx. 



(7) 



Let us calculate the spatial density n(x) for the poly- 
mer case. To do so we need to calculate the integral (|6]) 
substituting the energy <j5j> 



n(x) = — [ 

27Th J Z-Hx 



dp x 



(*)ex P (Ssin 2 (^))- ' 



(8) 



where Z{x) — exp /3 (ji — 
© can be replaced by a geometric series 



The integrand of 



n ( x ) = 7TT 
2irn 




Z(x) e 




d Px . (9) 



The sine function in the argument of the exponential can 
be written by using an identity, then we interchange the 
summation and the integration 



1 '^sj 



dp x - 



(10) 
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The last integral can be recognized as a modified Bessel 
function of first kind [60j , if we introduce a regulator for 
the integral that corresponds to the polymer length A 
15], with the result 



n(x) 



1 Twr 
2^h~ 



3 = 1 



V 4mA 2 



(11) 



We obtain the total number of particles by integrating 
(fTTj) over all space. The only x-dependent function is 
Z(x), and then the corresponding integration is straight- 
forward 



1 / 2tt ^ 



e 4ml' 



Io 1^7x2 7I7a ■ ( 12 ) 



A u V4mA 2 / jV2 



When A -C 1 , we can use the as ymp totic expansions of 
the modified Bessel functions [BO, Elf 



Io(u) 




(13) 



Consequently, the total number of particles is given ap- 
proximately by 



oo 



kT 

k=0 



-\ 2 ) k 



t (2mnT\ k r (fc + I) 



k y 

(14) 

where g v (z) are the so-called Bose-Einstein functions, 
being z = e 13 ^ the fugacity (58[ . Series (Tl4"|) is not conver- 
gent for arbitrary values of the polymer length. However, 
the approximation (fl~3|) is valid only for A <C 1. We no- 
tice that the leading term turns to be gi(z), the standard 
result [HI [59| . The second term can be considered as a 
correction of order A 2 , and is proportional to 172(2) 



N 



kT 

hw x 



gi(z) + x 



, mnT 



g 2 (z) + 0(X 4 ) 



(15) 



It is well know that a one-dimensional condensate can- 
not exist in the thermodynamic limit, due to the diver- 
gent behavior of the Bose-Einstein function 31 (1). Nev- 
ertheless, if we take into account the ground state energy 
of the system, we are able to obtain a well defined con- 
densation temperature [H, [H, [5{| . The ground state en- 
er gy a ssociated with our system was already calculated 
in [lOj, |l_5j, ll6| and is given by 



t d 2 / A 4 /4d 4 
e = ^-(l + _a [^ r 



(16) 



where d 2 — h/muj x , is the characteristic length of the 
one-dimensional oscillator, and a n (x) is the Mathieu 
characteristic function (60j . However, in the limit A -C 1 
one can regain the usual ground state energy plus correc- 
tions caused by the polymer scale 



huj x 



A 2 
32 



0(A 4 ). 



(17) 



At the condensation temperature, the value of chemical 
potential (i takes the minimum energy associated to the 
system. Substituting (p~7|) into (fT5]) and by using the 
properties of Bose-Einstein functions when eo / kT -C 1 , 
[58j we can reexpress (1151) at the condensation tempera- 
ture as follows 



N = 



hijj x 



gi (z c ) + X 



, mnT c 

~2h 2 ~ 



g2(z c ) + 0(X A ) 



(18) 



where z c = e -e °/ KTc . 

Setting A = we recover the usual expression for the 
number of particles [HI, 135] 



N= KT ° ln ( 2KT ° 



(19) 



where To is the usual condensation temperature in one 
dimension without polymer modifications. By using an 
iterating procedure 49], the condensation temperature 
turns to be 



kTq sa Hlu x 



N 



In 2 AT 



(20) 



Bose-Einstein condensates in one or two dimensions have 
been extensively studied [34| [3^, [59| (and references 
therein). The one dimensional condensation is seemingly 
not possible in the thermodynamic limit, in other words, 
the condensation temperature tends to zero when the 
number of particles tends to infinity [3J, l35l [54j, l55l . [59J . 
Nevertheless, if one takes into account the associated 
ground energy, the condensation is possible at finite tem- 
perature. Finite size effects arc needed in one dimen- 
sional systems, in order to make the condensation possi- 
ble. 

In order to calculate the condensation temperature T c , 
from (|18[) , we introduce the following well known expres- 
sions for the Bose-Einstein functions [58| 



ln(KT c /e ), 

C(2) - ^ 
£0 



1 +ln 



kTc 



(21) 
(22) 



After some algebraic manipulation, we are able to express 
the number of particles as follows 



N = 



1 - 



7 A 2 



1-C(2) 



A^ 
2d 2 
kT, 



7 huj q 



+ 0(A 4 ), (23) 



where A 2 = 2ttH 2 / 'mKT c , is the standard one dimensional 
thermal wave length. Notice that the leading term is 
identical to (fl9"|). with T c instead of T . From ([25]). we 
can notice that there are two kinds of corrections to the 
number of particles, associated with the polymer length. 
One is due to the ratio A 2 /d 2 , that is related to the effec- 
tive size of the ground state of the harmonic trap. The 
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second kind of correction is due to A 2 /A 2 , this term can 
be interpreted as a pseudo-interaction within the system 

To calculate the shift caused by the polymer scale A 2 , 
let us expand (|23l) at first order in A 2 = and T c = To, 
recalling that To is the usual condensation temperature 
given by equation ([2"0")l . The resulting shift AT C is 



T c -T 
To 



at; 

To 



A2_ 

2d 2 



\n2N-l 

o 



C(2)jV 
ln2JV 



1 + In 2N 



(24) 



Notice that the possibility to obtain a measurable co- 
rrection associated to the polymer scale A 2 (<5T C A ) re- 
quires that, if A(T C ) is the experimental error, then 
A(T C ) < |5T C A |, which in our case this entails 



A(T C ) 



< 



A 2 



2d 2 



In 27V 



C(2)jV 
In 2N 



1 + In 27V 



(25) 



In order to obtain representative bounds associated with 
the polymer scale, let us appeal to the current high pre- 
cision experiments for fgK in a 3-Dim condensate [31] . 
Although, from the experimental point of view there is no 
a real one dimensional condensate, it is possible to con- 
struct a quasi one dimensional condensate just by using 
extremely anisotropic traps 54]. Then, in principle, the 
use of this experimental accuracy is justified to bound the 
polymer scale A 2 . The shift in the condensation tempera- 
ture respect to the ideal result, caused by the interactions 
among the constituents of the gas is about 5 x 10~ 2 with 
a 1% of error 1311 ) . For instance, under typical condi- 



tions the number of particles is about N ~ 10 — 10 in 
three dimensions [34] : for one dimensional systems, the 
number of particles can be estimated up to TV ~ 10 4 
[52l [H, UK]. Using typical frequencies of order 100 Hz, 
allows us to bound the polymer scale up to A 2 < 10~ 16 
m 2 , which is notable. Indeed, as far as we know, this is 
the first bound associated to A, from low energy, earth 
based experiments. 

Notice that, if we increase the number of particles for a 
fixed trap frequency of order 10 2 Hz, the bound associated 
to the polymer scale decreases. On the other hand, for 
a fixed number of particles of order ~ 10 4 , again the 
bound decreases, when increasing uj. These facts suggest, 
that finite size systems are required in order to obtain 
representative bounds for A. 



IV. DISCUSSION 

Quantum gravity phenomenology suggests that the 
quantum structure of space-time could have some effect 



on the dynamics of particle motion. These effects can be 
modeled, for instance, as deformations of dispersion re- 
lations. Particularly, the quantization that arises in loop 
quantum gravity applied to the motion of a quantum par- 
ticle in one dimension (known as polymer quantization), 
gives a regularized Hamiltonian together with the intro- 
duction of a length parameter known as polymer scale. 
There are many proposals in the experimental search of 
these effects, however, is of special interest to study the 
Bose-Einstein condensates as they experimentally offer 
high accuracy measurements, which allows to establish 
bounds to the polymer length scale. 

The main goal of this paper was to constrain the poly- 
mer length parameter by studying an non-interacting 
one dimensional polymer Bose-Einstein condensate in an 
harmonic trap within the semiclassical approximation. 
We have proved that the corresponding expression (Ti~4"[) 
generalize the so-called Bose-Einstein functions that ap- 
pears in the standard case, to a series in terms of powers 
of the polymer length scale. Using typical experimental 
values for the number of particles in the condensate, the 
mass, and the frequency of the harmonic trap, we were 
able to establish a bound for the polymer length up to 
A 2 < 10~ 16 m 2 . To improve this bound, finite size sys- 
tems in one dimension are required. 

Notice that if we consider the ratio (A/0?) 2 as an ef- 
fective dimensionless parameter then in principle, we are 
able to compare with the deformation parameters sug- 
gested in other approaches [23T - [25j by using condensates. 
We can see that the order of magnitude are quite similar 
and opens the possibility to relate different scenarios as 
was recently suggested [57] . 

It is possible to consider more general traps like power 
law potentials, the 3D case, and interactions among the 
constituents of the gas, to improve further the bound on 
the polymer scale, and will be presented elsewhere [26] . 
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